Interaction effects on ID fermionic symmetry protected topological phases 
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In free fermion systems with given symmetry and dimension, the possible topological phases 
are labeled by elements of only three types of Abelian groups, Zi, Z2, or Z. For example non- 
interacting ID fermionic superconducting phases with Sz spin rotation and time-reversal symmetries 
are classified by Z. We show that with weak interactions, this classification reduces to Z4. Using 
group cohomology, one can additionally show that there are only four distinct phases for such ID 
superconductors even with strong interactions. Comparing their projective representations, we find 
all these four symmetry protected topological phases can be realized with free fermions. Further, we 
show that ID fermionic superconducting phases with Zn discrete Sz spin rotation and time-reversal 
symmetries are classified by Z4 when n = even and Z2 when n — odd; again, all these strongly 
interacting topological phases can be realized by non-interacting fermions. Our approach can be 
applied to systems with other symmetries to see which ID topological phases can be realized by free 
fermions. 



Symmetry protected topological (SPT) phasesfTJE] are 
short range entangled states with symmetry protected 
gapless edge excitations. [3HH] The Haldane phase on a 
spin-1 chain p, [TU] and 2D/3D topological insulators [TT]- 
[TS] are examples of SPT states. Using K-theory or 
topological terms, all free-fermion SPT phases can be 
classified [T71 [T5] for all 10 Altland-Zirnbauer symmetry 
classes of single-body Hamiltonians. It turns out that 
different free fermion SPT phases are described by only 
three types of Abelian groups, Zi, Z2, or Z. 

With interactions the classification is more varied, 
however we must first describe the symmetry differently. 
Instead of specifying the symmetry of single-body Hamil- 
tonians, we treat the free fermion systems as many-body 
systems and specify the many-body symmetry of their 
many-body Hamiltonians. Only in this case can we accu- 
rately add interaction terms to the many-body Hamilto- 
nian that preserve the many-body symmetry, and study 
their effect on the SPT phases of free fermions. A classi- 
fication of various free-fermion gapped phases given their 
many-body symmetry can be found in Ref. 1201 

Fidkowski and Kitaev (also Turner, Pollmann and 
Berg) studied interaction effects in one case: In 
their ID time-reversal (TR) invariant topological 
superconductor, [5TH13] the Z classification in the free 
case breaks down to Zg with interactions that preserve 
TR symmetry. Here we present another model which il- 
lustrates the effects of interactions on ID SPT phases. 
We begin with a lattice Hamiltonian for a ID supercon- 
ductor with both TR and spin-rotation symmetries, 
described by the Z classification in the free case. With 
the addition of weak interactions that preserve these sym- 
metries, the classification reduces to Z4 (see Table |l]). 

We compare these four fermionic phases to the four 
phases predicted separately from group cohomology |i4:| 
m [21] (a method valid for strong interactions). We 
find each fermionic phase has a distinct projective 
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TABLE I. Symmetry groups described by the ID Hamiltonian 
in Eq. 1 (where Z2 is generated by the time reversal transfor- 
mation), with their free fermion classification and how they 
reduce with interactions. The latter remains true with strong 
interactions so all phases can be realized with free fermions. 



representation^ i4j and since group cohomology also 
gives rise to four and only four distinct phases, [7] we 
conclude that all free fermions can realize all strongly 
interacting SPT phases in this case. We further study in- 
teraction effects on a ID superconductor with Z„ discrete 
Sz spin rotation and TR symmetries. For this symmetry 
group, we find the SPT phases are classified by Z4 when 
n = even and Z2 when n = odd. Again, these results are 
separately obtained both from perturbing our fermionic 
lattice Hamiltonian and from the group cohomology clas- 
sification for strong interactions - showing that again, all 
the strongly interacting topological phases can be real- 
ized by non-interacting fermions. 

Free fermion lattice model — We write a ID Hamilto- 
nian with a trivial and two non-trivial phases 



H 



(1) 



where the first term is typical nearest-neighbor hopping, 
the second term Ag represents on-site pairing and the 
last term with Ap pairs electrons on adjacent sites. 
This Hamiltonian satisfies time-reversal T and Sz spin- 
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FIG. 1. Phase diagram when varying the parameters and 
Ap-. The phase boundaries are As = ±Ap which separate 
three phases denoted by — +1,0 and -1. In the center phase 
we can make As arbitrarily large without closing the gap. 
That limit describes on-site pairing where any cut through 
the system separates both parts cleanly without leaving edge 
states - allowing identification of this trivial N = phase. 



rotation symmetries specified on cf^ = {ci^, 64} as 



-ie/2 



ie/2 



so that THf-^ = H and 



-He 



-ies. 



H. As the 



bandgap closes to leave just the hopping component when 
As = ±Ap, we obtain the phase diagram in Fig. [T] 

We start by identifying the trivial phase = 0: When 
I As I > |Ap|, we can arbitrarily increase the strength of 
As without closing the gap. In the limit of As much 
larger than the other terms, the Hamiltonian simply re- 
duces to on-site pairing — any cut through the system 
separates both parts cleanly leaving no boundary states 
(see Fig. [T]): this is the trivial = phase. 

Next, we look for ground-state degeneracy at the inter- 
face between this phase and its neighbors. This is most 
conveniently done in a low-energy continuum model, 
where the effective Hamiltonian becomes 



H = -i dx¥ 



in a basis of right and left-moving fermion operators 
= (V'fl.-t-, iV'li'^j iV'l^, V'Lt) close to the Fermi surface. 
Here m = Apl — A^az- 

Smoothly varying our mass term rn(a;) across an in- 
terface, we set Ap(x) = ^(1 + tanhx) and As(x) — 
|(1 — tanhx). This has the zero-energy solution 



(2) 



This complex fermion operator ("00+ 7^ ^0+ ) with energy 
E = contains a double degeneracy (empty or filled) that 
allows labelling of Ap > |As| as the non-trivial A^ = 1 



phase. This mode transforms under symmetry as 




(3) 



Since the two degenerate states differ by S'^ — 1/2 and 
are related by time-reversal, each state carries quantum 
number of Sz = ±1/4 respectively. 

Using the symmetry relations in Eq. [3] we check if any 
perturbations in the Hamiltonian can shift the energy of 
this mode. We find density terms SH — c'0j_|_'0o+ are 
forbidden by TR, hence our ground-state degeneracy is 
protected by system symmetries — this A^ = 1 phase is 
stable against perturbations. 

To find the A^ —1 phase, we change Ap(a;) -> 
—Ap{x) and upon repeating our procedure, find a dif- 
ferent zero mode solution that we label 



■00- 



dxsech(x)(i^^ 



(4) 



Cicr and instead transforms as 



T 



■00- \ rf^-l „ { 00- 



-ie/2 



■00- 



(5) 



This state has stable ground-state degeneracy as 5H = 
cV'o-V'o- is also forbidden by TR, indicating Ap < |As| 
is a non-trivial phase as well. 

Is it meaningful to label this second non-trivial phase 
A^ = — 1? We examine what happens upon stacking two 
chains both with non-trivial phases but the first with 
Ap > |As| and the second with Ap < |As|. Thus the first 
chain would have a zero mode ^0+ and the second would 
have V'o-. We find the coupling 5H = c'ip\_^ij)a- + h.c. is 
allowed within system symmetries and makes the ground 
state nondegenerate. So two chains with the two dis- 
tinct zero modes (labelled -f and — ) combine to become 
trivial, indicating the two phases should be labelled with 
opposite index. Naturally then the phase with i/^q- would 
be the A^ = — 1 phase, so this model indeed gives three 
symmetry protected phases A^ = —1,0 and -t-1. 

While two chains containing zero modes with opposite 
index become trivial, we further consider the stability of 
two chains containing zero modes with the same positive 
(or negative) index. This may generalize to larger inte- 
gers in the Z group, so we now examine the stacking of 
two chains with similar index more systematically. 

A generic coupling term (see Fig. [2]) is 6H = 
cip!^g^Mab4'+b + li-C. Here a and b are indices running 
over the chain number 1,2... e.g. ip^i denotes a zero 
mode from the A^ = 1 phase in the first chain; and Mat 
is any generic coupling between these two operators. We 
examine the simplest case of a = 1 and 6 = 2. Terms 
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FIG. 2. We stack two chains in the same non-trivial phase 
with positive (or negative) index to see if their edge states are 
stable. With the first chain a = 1 and the second 6 = 2, Mat 
is any coupling between them. We find all possible couplings 
are forbidden by our system symmetries so two similar modes 
are stable and form the N = 2 phase. 



FIG. 3. Without interactions V1234 = the ground-state 
degeneracy for four zero modes is 2^ = 16. With interactions 
V1234 / two states are split by 5E — ±|Vi234|, making the 
ground state nondegenerate and this A'^ = 4 phase trivial. As 
= 4 is now smoothly connected to the trivial A'^ = phase, 
our classification reduces from Z to Z4 with interactions. 



of the form 5H = cV'^]^Mi2'0+2 + h.c. are forbidden by 
TR symmetry as specified in Eq. |3j while fermion pair- 
ing terms such as 5H — c'!/'^j^Mi2V'+2 violate 
spin rotation symmetry. As there are no other quadratic 
fermion terms, the stacking of two chains is stable against 
perturbations and combine to give an A'^ = 2 phase. 

Hence adding a number of ID chains with positive in- 
dex corresponds to a positive integer in the Z group. The 
negative numbers are obtained simply by stacking chains 
with V'o- • As we showed earlier that a pair of ■00+ ^nd 
'0o_ coupled together become trivial, the integer N in 
our Z group is the difference between all positive and 
negative zero modes. Then for each phase labelled by A^, 
the ground-state degeneracy is 2l^l. 

Interaction effects — Now we allow couplings with an 
arbitrary number of fermion operators. We look at terms 
with four and two operators which take the general form 

5H = VahcdV'+a^+&V'+cV'+d + W^abV'+a'/'+b + ^.C. (6) 

5H is compatible with both TR and Sz spin-rotation 
symmetry when Vabcd and Wab satisfy certain conditions. 

A possible term couples four separate chains through 
an interaction with only V1234 ^ 0. Such an 5H is in- 
variant under both TR and spin-rotation symmetry. 
This couples the two states |0101) and |1010) in our four- 
mode basis (0 and 1 denote unoccupied and occupied 
respectively for each of the four chains). Without inter- 
actions we have a ground-state degeneracy of 2"^ = 16; 
with interactions two of these 16 states split in energy by 
5E = ±1^1234!, see Fig. [s] This makes the ground-state 
nondegenerate and the phase A^ = 4 trivial. 

Since four chains with all positive (or negative) index 
are equivalent to the trivial phase, we can smoothly con- 
nect the A^ = 3 phase to the A^ = — 1 phase by adding 
four chains with all negative index. So with only three 
distinct non-trivial phases, the Z integer classification for 
free fermions reduces to Z4 in the presence of interactions. 

Four- fermion interaction terms also reduce the ground- 
state degeneracy in the N = 2 phase from 2^ = 4 to a 



two-fold degeneracy. The term 

SH = V^1122('0+i^/'+l - ^+l1pli){lpl_2'^+2 - '0+2V'+2) 

causes two states |00) and |11) to shift in energy by V1122 
while two other states |01) and |10) shift by —^1122- As 
we still have doubled ground-state degeneracy, the state 
N = 2 remains non-trivial. To summarize, interaction 
effects reduce our degeneracy leaving three non-trivial 
phases each with a two-fold ground state degeneracy. 

Four distinct projective representations — Our results 
depend on perturbating from free fermions, i.e. we check 
the stability of free fermion phases with weak interac- 
tions. This method may not capture all possible inter- 
acting phases as strongly interacting topological phases 
may not adiabatically connect to free fermion phases. Or, 
different phases from weak interactions may become the 
same phase with strong interactions. To address these is- 
sues, we calculate the projective representation for each 
phase to see if they are truly distinct. 

Using the symmetry operations defined in Eqs. [3] and 
[5] we write their matrix representation on the degenerate 
subspace. In the A^ = 1 phase with basis |0+) and |1+) 

Ug ^ M{Ub) = ^ ^_,^g/2 ^ , M{f)K = a,K 

where T = U-j^T, a rotated TR operator we can intro- 
duce since Ug and T commute. Then 

M{f)KM{Ug) = e'^^^M{Ug)M{f)K (7) 

and M{f)KM{f)K = 1. This is a projective repre- 
sentation as the phase in Eq. [7] cannot be removed by 
adding any phase factor to M{Ug) that changes the spin 
quantum number. 

Moving to the A^ = — 1 phase with ground states |0_) 
and we have the same representation for M{Ug) 

while M {f)K = e''-''^'' I'^u^K in this case. Eq. [rjremains 
true but now M {f)KM{f)K = -1 so we have a different 
projective representation. 

In the N = 2 phase, our ground states are four-fold: 



4 



|0+0+),|0+l+),|l+0+) and |1+1+). Here 

We) - j « ^ j , 

M(f )if = (7,^(^ ayK. (8) 

While M{f)K and Af([/e) commute this time, 
M{T)KM{T)K = — 1 again making a third non-trivial 
projective representation. 

As each non-trivial phase has a distinct non-trivial pro- 
jective representation, they remain distinct phases even 
when interactions are strong. We can compare our results 
to the unperturbative bosonic classification in ID ob- 
tained by group cohoniologv[^ 171 as we can bosonize 
our fermionic model. The resulting bosonic model would 
have the same symmetry U{1) x Zj with phases classi- 
fied by Z2 X Z2, i.e. four distinct projective represen- 
tations there correspond to four different strongly inter- 
acting phases. Our fermionic results similarly contain all 
four phases with these distinct projective representations, 
so this model realizes all possible non-trivial phases with 
strong interactions. 

Modifying symmetry from U{1) to Z„ spin-rotation — 
As our fermion model respects Sz spin-rotation and TR 
symmetry, it naturally contains Z„ discrete spin-rotation 
as well. We can replace U{1) spin-rotation by Z„ spin- 
rotation, i.e. rotation by an arbitrary angle is now con- 
strained to values oi 6 — li^jn and our new symmetry 
group has generators time-reversal T and discrete 6*2 ro- 
tation R= e}^'^ satisfying 

y2 ^ ^« ^ {~)^'' . TR = RT. (9) 

Here (— is the fermion number parity operator. 
When n — even, this group G{T, Zn) is generated by 
R and f = K^/'^T, so G(r, Z„) = Z2„ x Zj. When n = 
odd, we find that R = RT alone generates this group 

G(T,Z^) = Zl„m- 

For n > 2, no new fermion bilinear terms are allowed 
so the free fermion classification does not change from Z. 
In the case of n = 1, new quadratic terms of the form 
5H = C'i/'^_]^^^2 + -ti-C- are permitted. This term couples 
two chains forming the N = 2 phase to make the ground- 
state nondegenerate. The N = 2 phase becomes trivial 
and the classification for n — 1 reduces to Z2. 

Similarly for higher n, we can always add interacting 
terms with 2n tp^ operators similar to the term in the 
n = 1 case above. For n = 2 for instance, this term is 
SH = c^+i'!/'+2''/'+3V'+4 + h.c. Such interactions couple 
2n zero modes each in the iV = 1 phase to render the 
ground-state nondegenerate. In effect, Z„ spin-rotation 
symmetry allows interactions that reduce the classifica- 
tion to Z2n- 

We had established that under U{1) spin-rotation sym- 
metry, interactions reduce the classification to Z4. In- 
cluding more interactions as allowed by Z„ spin-rotation 



further reduce the classification to Z2„. Taken together, 
we find there is no effect on even n which remains Z4 
since 2n is a multiple of 4. Odd n, however, reduces 
to a Z2 classification (as 2 becomes the largest common 
denominator between 2n and 4). 

The number of non-trivial phases can be compared 
to and matches with the group cohomology prediction 
Z2 X Z2 for even n and Z2 for odd n[5S]. We find that 
different symmetry groups with the same free fermion 
classification reduce to various results (here Z4 or Z2 
are examples) in the presence of interactions (summary 
in Table |l|. As verified by comparison of these phases 
with group cohomology, all possible strongly-interacting 
phases can be realized by free fermions in this model. 

Lastly, we note that our classification is protected only 
by system symmetries of spin-rotation and TR. As shown 
earlier, without such symmetry a term SH = ctpl_^_'!pQ^ 
would be permitted which renders the ground-state non- 
degenerate and the classification trivial (Zi). 

Discussion — We study the SPT phases of ID 
fermionic superconductors with TR and Sz spin-rotation 
symmetries. If fermions do not interact, their classifica- 
tion is given by the Z group; with weak interactions this 
reduces to a Z4 classification. As each of our four fermion 
phases have distinct projective representations, they cor- 
respond to four distinct phases by comparison with group 
cohomology, which predicts four and only four different 
gapped phases even with strong interactions. Hence all 
distinct symmetric gapped phases with strong interac- 
tions are realized by non-interacting fermions in this case. 

As the edge states in our ID superconductors are de- 
scribed by complex fermions which have double the de- 
grees of freedom of Majorana fermions, it is unsurprising 
that the interacting classification is half of the result from 
Kitaev and Fidkowski's Majorana modei pm23j . 

We further studied the SPT phases of ID supercon- 
ductors with TR and Z„ discrete Sz spin-rotation sym- 
metries, to find they are classified by Z4 when n = even 
and Z2 when n = odd. Again, as phases in our fermionic 
model matches with the group cohomology prediction, 
all gapped phases of these ID fermionic superconductors 
are also realized by non-interacting fermions. 

Interactions on different symmetry groups with the 
same free fermion classification give rise to varied re- 
sults (Table |l|. Here perturbing from a free fermion 
model gives all strongly interacting phases, however in 
other cases such phases may not be realizable with free 
fermions. Lastly, the effectiveness of this method remains 
open especially in higher dimensions where additional 
tools may be needed. Further study of different symme- 
try groups or in higher dimensions would be worthwhile. 

Towards the completion of this paper, we noted the 
work of A. Rosch (arXiv:1203.5541 ) which shows "a topo- 
logical insulator made of four chains of superconducting 
spinless fermions characterized by four Majorana edge 
states can adiabatically be deformed into a trivial band 
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insulator" via "interactions to spinful fermions", which 
has some relation to our Z4 classification of ID fermionic 
superconducting phases with TR and Sz spin-rotation 
symmetries. 

We thank Zheng-Xin Liu, Andrew Potter and Xie Chen 
for helpful discussions. This work is supported by NSF 
Grant No. DMR-1005541 and NSFC 11074140. 



Free fermion classification 

Given our many-body symmetries T and Q (where Q 
denotes our U(l) symmetry), the mass matrix M should 
satisfy the following relations [20 



{7,T} = {7,Af} = {r,Af} 
[g,7] = [Q,T]=0 



0; 



1 



(10) 



where 7 is simply the hopping part of our Hamiltonian. 

It is simplest to choose a basis where Q — I®e (here we 
write £ — —iCy) as the the anti-symmetric Hamiltonian 
we call A takes the form [501 

A^Ha®I + Hs®e; H = Hs+iHa (11) 

and Hg — Hf and Ha = —H^; i.e. the space simplifies 
into the complex space Cq. 

Given this choice our operators can take the form 



7 = cr^ 



K)/(8)/ 
> /(g) e 
g/(g)e 
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(12) 

Our Hamilto- 



which satisfy the conditions in Eq, 
nian becomes the expression given in Eq. [TT] so we easily 
identify the space of our mass matrix: 



M = £ (g) {/,n • cr, -/} (g / 



(13) 



i.e. the three possible phases have eigenvalues of -f 1, ±1 
and —1 respectively. (Note that here n a = nxcrx + nye + 
Uzdz i.e. M is real.) These three phases correspond to 
the = — 1, and 1 phases discussed in the paper. 

Note the fermion number operator N — J^ii^l^^it + 
cl^Cii) in this basis is A = g) cr^ g) £. As this does not 
commute with the mass matrix in both the ±1 phases, 
fermion number is not conserved — for our basis we 
choose U{1) symmetry to represent spin conservation in- 
stead of particle number conservation. 



Group cohomology classification 

Here we calculate the group cohomology for the sym- 
metry group we denote as G{T,Zn), generated by the 
time reversal T and discrete Sz rotation R = e''^^^. 
The generators T and R satisfy 



= (-) 



Nf 



R- = (-) 



Here (— )^^ is the fermion number parity operator. 

When n = even, we see that G{T, Z„ ) is generated by 
i? and f = R'^^^T, which satisfy 

f 2 = 1, i?" = (-)^^, TR = RT. (15) 

Thus G(T, Zn) = Z2n X when n — even. This allows 
us to obtain [7] 

H^[G{T,Zn),UTil)]=ll, n = even. (16) 

When n = odd, let R = RT. Wc find that R" = 
[(_)A^F]i+Vy and -R^n = (-)^^. So a combination 
of i?" and i?^" will give us T and a combination of T 
and R will give us R. Thus R alone generates G(T, Zn) 
when n = odd. We find that, for odd n, G{T, Z„) = Zj^^ 
which is generated by R with i?*" = 1. In Zf^, R^ is 
anti-unitary if fc = odd, and unitary ii k ^ even. We 
find H'^[ZI^,Ut{1)] = Z2 (at least for n = 1,2,3,4,5,6). 
Thus 



F2[G(r,Z„),C/T(l)] =Z2 



odd. 



(17) 



Thus the symmetry protected topological phases with 
time reversal and Zn discrete 5*2 spin rotation symme- 
tries are classified by for even n and Z2 for odd n. 



TR = RT. 



(14) 
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